Pull-In Range of a Phase-Locked Loop 
With a Binary Phase Comparator 

By JAMES F. OBERST 

(Manuscript received May 25, 1970) 

We develop a method for calculating the pull-in range of a phase-locked 
loop with a binary phase comparator and an arbitrary loop filter. Complete 
numerical 7-esults are presented for loop filters of the phase-lag and low-pass 
types. The problem of stability is also considered, and it is proved that with 
these loop filters no steady-stale phase jitter can exist after frequency 
acquistion has been achieved. 

I. INTRODUCTION 

The phase-locked loop (PLL) is an important element of many- 
modern communication and control systems. A PLL block diagram is 
shown in Fig. 1. The input i>i(<d £ + #i) is a narrow-band signal with 
carrier frequency <o and phase 6\{t). This phase is compared with 
the phase 6 2 (t) of the voltage-controlled oscillator (VCO) in the phase 
comparator (PC). The PC output /(<£), where = $ x — 0*, is filtered 
by the loop filter H{p) and applied to the VCO control terminal. 

Depending on the values of the PLL parameters, the phase error <f> 
can be kept small even with input phase modulation. Thus with 6\{t) 
= Clt + 0io, which represents a constant input frequency offset, the 
system can produce a synchronized signal v z (t) with frequency w + 
ft. This synchronization capability leads to PLL applications in carrier 
extraction, 1 frequency synthesis,- narrow-band filtering, 3 FM demodu- 
lation, 3 timing extraction in PCM and data transmission systems, 4 etc. 

In this paper, we examine the acquisition, or pull-in, range of a 
PLL with a binary PC. We present numerical results for the special 
case of a second-order PLL with either low-pass or phase-lag loop 
filter. 

The PC characteristic considered here is the binary curve shown in 
Fig. 2. It is of interest in at least three situations. First, since many 
synchronization systems are designed to operate with very small 
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Fig. 1 — PLL block diagram. 

phase errors, dynamic range limitations in the PC circuitry often pro- 
duce severe saturation effects. The characteristic of Fig. 2 corresponds 
to the extreme case of vanishing linear range near zero phase error. 
However, it is a useful approximation for systems with small but non- 
zero dynamic range for the purpose of studying pull-in performance. 
Second, a binary PC can be easily implemented with logic circuits. 
The resulting characteristic differs from the ideal of Fig. 2 by exhibit- 
ing small hysteresis about the zeros at <£ = titt, but this hysteresis has 
no effect on the pull-in range achieved. Finally, J. J. Stiffler 5 has 
shown that for a first-order PLL with additive white gaussian noise 
and no frequency offset, the cross-correlation type PC which mini- 
mizes Pr { | | > <£o } for all <f> has the characteristic of Fig. 2. 
Although such a square-wave correlation function is unrealizable, 
this result suggests that PLLs employing other types of PC having 
this characteristic are worthy of consideration. In addition, similar 
"bang-bang" control characteristics are known to be optimum for 
PLL acquisition. 6 



II. CALCULATION OF PULL-IN RANGE 

The phase model corresponding to the PLL block diagram in Fig. 
1 is shown in Fig. 3. We assume that the gain of the loop filter H(p) 
is unity at DC. The input-signal frequency differs from the VCO center 
frequency by O rad/s: 

6,(1) = Qt + 1O . (1) 

It is convenient to normalize the detuning n to the dc loop gain* <*: 

7 = Q/«. (2) 



* Since no gain can be defined for the binary PC being considered, the symbol 
a does not represent the usual small-signal loop gain. 
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Fig. 2 — Phase comparator characteristic. 

When the normalized detuning y is not too great, the VCO frequency 
changes toward the input frequency, and eventually the PLL synchro- 
nizes to the input signal with zero frequency error and finite phase 
error. The normalized lock range y L is the maximum detuning | y \ 
for which the PLL can remain locked after synchronization has been 
established. Inspection of dc conditions in the phase model of Fig. 
3 shows that y L = 1. The normalized pull-in range y P is the maximum 
| y | for which eventual synchronization is assured from any initial 
conditions in the loop filter and VCO. In general y P < 1 for PLLs of 
order higher than first. The order of a PLL is defined as one plus the 
number of poles in H(p). Calculation of y P is the subject of this paper. 
The method employed here is similar to that used by A. J. Gold- 
stein 7 to calculate y P for a PLL with a sawtooth PC. Due to the binary 
nature of /(<£), the PC output waveform f[<f>(t)] is piecewise constant, 
assuming only the values ±1. Assume that the PLL is not synchronized 
to the input signal. Then <f>(t) increases with time (for n > 0), and the 
waveforms <f>(t) and f[<f>{t)] appear as shown in Fig. 4. The time ori- 
gin has been selected so that <f>(0) = 0. The transition instants are 
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Fig. 3— Phase model of PLL. 
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Fig. 4— Phase and PC output waveforms. 



where 

</>('» i) = ( 2 J _ IV (negative transition), 

ff,(t j2 ) = 2 pc (positive transition). (4) 

An expression can be written for /[0(f)] by summing all of the seg- 
ments corresponding to 2?r increments in <f>(t): 

/[*(<)] = E [«(' - */-!.■) - 2w (* ~ '») + w (' ~ ^)]< (5) 

| =-oo 

The ;' = 1 segment is shown crosshatched in Fig. 4. Since the PLL is 
not synchronized to the input signal, the steady-state PC waveform 
f 8 a[<t>(t)] is periodic, and the transition instants can be written: 



(6) 



t ia = j(T 3 + T t ), 

tn = KT 3 + T 4 ) - T 4 . 

T a and T 4 are the times between transitions as indicated in Fig. 4. 
Using equation (6), equation (5) becomes: 

/..[*(<)] = E W - [j - i][r, + rj) 

1- -oo 

- 2u(i - j[T 3 + rj + T 4 ) + u(* - ;[T, + rj)]. (7) 
From equation (1) and Fig. 3, the phase error in steady state is: 



0..(O = Qt + 0o - /..[*(0] * £ _1 {^} . 



(8) 



where 4> is some constant. Since } 88 [<f>(t)] is composed only of step 
functions, the last term in equation (8) can be written as a sum of 
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delayed functions g(t), where 



JBlrfOI -W-^P- W 



The general expression for SS (£) is then 

0„(O = m + «o - a E fo(* - [j - !][ r 3 + rj) 

j •» — CO 

- 2 9 (t - j[t 3 + rj + r 4 ) + 0(f - icr. + rj>]. OT 

From equation (4) and Fig. 4, we have the following conditions on 

*..(0) = 0, 

*..(r 8 ) = tt, (ii) 

*..(Ts + ^ - 27T. 

These conditions are sufficient to determine the unknown constants <f> , 
T a , and T 4 in equation (10). Thus equations (10) and (11) together 
define the relationship between the normalized detuning y and the 
loop-filter parameters [through g{t)] which must be satisfied for the 
PLL to be out-of-lock in the steady state. Then clearly y P is the 
minimum value of y for which these equations possess a solution. 

III. RESULTS FOR SECOND-ORDER PLL 

In this section, the method derived above is applied to the second- 
order PLL with loop filter 

H(p) _ L±^£. (12) 

This is a phase-lag filter for < T 2 < Ti . It becomes a simple low- 
pass filter for T 2 = 0, and setting T 2 = T x reduces the PLL to first 
order. The corresponding g(t), from equation (9), is 

g(t) = [t- CZ\ - T 2 )(l - exp (-t/T^Mt). (13) 

In Appendix A, equation (10) is rewritten using equation (13) and is 
evaluated at t = 0, T 3 , and T a + T 4 . Equation (11) is then applied, 
along with the normalization 

r, - aT i} i = 1, 2, 3, 4. (14) 
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The equations which result are 

2tt + t 3 — T 4 
7 = r, + r. ' 

4(T, - T 2 )(t 3 + T 4 ) 



(15) 



= [r 4 (r 3 + w) + r 3 (r 4 - x)] [coth £■ + COth ^J ■ (16) 

According to the discussion following equation (11) , the pull-in range is 

27T + T 3 J— T4 r.-v 

7 P = min — , (17; 

Ti.T«>0 T 3 T r 4 

subject to the constraint equation (16). It is important to note that 
equation (15) is simply the dc balance equation for the PLL model 
of Fig. 3, and therefore holds for any H(p) with unity dc gain. Since 
equation (15) gives y explicitly in terms of t 3 and t 4 , it can always 
be used to eliminate y from a constraint equation corresponding to 
equation (16). Therefore y p can be calculated from equation (17) 
for any H(p) subject to the appropriate constraint equation which 
relates the loop-filter parameters to the transition-time parameters 
t 3 and t 4 . Hence only <£ 8S (0) and <fe„(T a ) actually had to be evaluated 
in Appendix A. 

The method employed to evaluate y p for various values of ti and t 2 
was to choose t 3 > 0, use equation (16) to obtain the corresponding t 4 , 
and calculate y from equation (15). The t 3 , t 4 relationship was found 
to be single-valued for all loop filters investigated. Examples of the 
behavior of y with t 3 are given in Fig. 5. The filter parameter r is 
denned as 

r = T 2 /T,. (18) 

In all cases the curves y(r 3 ) are smooth and exhibit a single local min- 
imum. This minimum is found by computing a sequence y(T 8n ), where 
T3n > t 3 ,„_i . When this sequence begins to increase, the minimum y p 
has just been passed, and can be estimated accurately from the last 
three computed values of y. 

Curves of y P versus t x with rasa parameter are presented in Fig. 6. 
Several characteristics are notable. First, y p = 1 for r ^ 0.5 indepen- 
dent of Ti . Second, as tj increases with r constant, y v approaches an 
asymptotic value y„ a which is a function only of r. Later we shall derive 
an explicit formula for y Pa . The same results are presented in a differ- 
ent way in Fig. 7, which shows curves of constant y v on the t, , t 2 
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Fig. 5 — 7 versus t 3 . 
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Fig. 7 — r 2 versus n with parameter y P . 

parameter plane. Below we derive the equation for the curve in this 
figure corresponding to y p = 1 . 

IV. FURTHER RESULTS 

Let us consider the important case of very large n, which corre- 
sponds to strong filtering in the PLL. Noting that 



Lim coth x = - 



and using equation (18) , equation (16) becomes for large n : 

4r,(l - r)(r 8 + u) = [r 4 (r 3 + x) + r 3 (u ~ *)] ^ + ^T I 
This simplifies to 



(19) 



(20) 



u = ^V- (2D 

2rr 3 + x 

Substituting equation (21) into equation (15) gives the result for y 
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with ti large: 



7 = 



Tt\ + 2rTTa + 7T* 
rT : , + 7rr :t 



All minima of y must satisfy 






(22) 



(23) 



Applying equation (23) to equation (22) yields a single positive value 
of t 3 : 



t 3 = 



-[■ - e 



- - 1 

r 



r < 0.5. 



(24) 



(25) 



1 - 2r 

Substituting equation (24) into equation (22) gives the result 

(2[r(l - r)]*, ^ r < 0.5, 

U, r £ 0.5. 

This agrees with the numerical results in Fig. 6 and with a result 
obtained by M. V. Kapranov by a different method in an untranslated 
Russian paper. 8 The existence of only a single minimum of y(r 3 ) for 
large tj supports our hypothesis of a single minimum for all ti which is 
based on the curves in Fig. 5. 

The region y p = 1 in Fig. 7 corresponds to r, , r 2 such that y ^ 1 
for all t 3 , t a . From equation (15), this implies that 

T 4 ^ 7T. (26) 

Thus from equation (16), t, , t 2 , and t ;h on the y p = 1 boundary must 
satisfy 



4(r, - r 2 )(r 3 + tt) = 7t(t 3 + jr)[coth^- + COth ^" ] • 



(27) 



Examination of Fig. 5 shows that the critical y curves approach y = 1 
from above for very large r :i . Using 



Lim coth x = 1 , 



equation (27) becomes 

4(t, — r,,)r ; , = 7T7- :t 



1 + coth — 
2r u 



(28) 



(29) 
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This reduces to the simple expression for the y p = 1 curve : 

r 8 - n - |/[1 - exp (-*/*)]. (30) 

Finally, let us consider the following problem. The periodic behavior 
of <f> ga (t) assumed throughout this paper which led to equations (15) 
and (16) is known as a limit cycle of the second kind in the phase 
plane. The nonexistence of such limit cycles for |y| < y p proves that 
frequency lock is eventually attained for these values of y. Physically 
speaking, <f> aa (t) cannot increase monotonically with time as assumed 
in Fig. 4. However, proper synchronization of the PLL requires that 
phase lock also be achieved. This means that after a long enough time, 
the system comes to rest: 

Lim 0(0 = 2ri7r, | 7 I < T P • (31) 

l-»00 

Because of the gross nonlinearity /(<£) in the system considered here, 
it is not obvious that equation (31) will be satisfied. Specifically, it is 
conceivable that a series of self-sustaining overshoots in tf>(t) could 
become established after pull-in which would produce a periodic phase 
jitter. This behavior corresponds to a limit cycle of the first kind in the 
phase plane. Although this problem is not solved in general here, a 
test which is valid for any H (p) is applied to the phase-lag filter case 
in Appendix B. It is found that in this case, phase lock described by 
equation (31) is always achieved. 

v. CONCLUSION 

A method has been presented for calculating the pull-in range y p of 
a PLL with a binary phase comparator and an arbitrary loop filter. 
The result is obtained as the minimum value of a function of two 
variables, subject to a constraint equation which relates these variables 
to the parameters of the loop filter. Complete numerical results for 
y p were obtained for loop filters of the phase-lag and low-pass types. 
Explicit formulas were given in this case for the asymptotic value of 
y p with strong loop filtering, and for the set of filter parameters which 
result in unity pull-in range. Finally, it was proved that no steady- 
state phase jitter can exist after pull-in with these loop filters. 
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APPENDIX A 

Evaluation of fa t (t) 

Consider one period of <f>„(t), and define the functions fa(t) as: 

♦..« = ><*• o***r» (32 ) 

k(0, t 3 ^ t $ r, + t 4 . 

From Fig. 4, 0i(£) includes all terms in equation (10) which correspond 
to transitions prior to t n = TV Using equation (13) in (10), <f>i(t) 
becomes: 

Mt) = at + O - a £ [i - (; - i)(T 3 + rj - (7\ - ro 

,--00 

•(1 - exp [-(I - y - 1][T 3 + rjyrj)] 
+ 2« £ [i - i(r 3 + r 4 ) + r 4 - or, - ro 

J --00 

• (i - exp [-(t - j[t 3 + rj + ro/rj)] 
-«EP- i(^a + ro - (2*1 - r a ) 

»--oo 

•(1 - exp [-(* - j[T 3 + rj)/r,])]. (33) 

Absorbing all constant terms into <f> , equation (33) becomes: 

<h(0 = m + <t>' Q - at 

- aCA - T 2 ) exp (-J/T,) £ [exp [j(T 3 + rj/TJ 

l—oo 

- 2 exp [(j[t 3 + rj - r 4 )/rj + exp y(r, + jy/rj]. 

= Q£ + 0o" - at - 2a(T, - T 2 ) 

. exD / ,/« ) 1 ~ exp (-TM 

exp( «/^) i _ exp [-(T 8 + TO/r,] (34) 

0o is eliminated from equation (34) using 

0i(O) = 0, (11a) 

which yields 

0,(0 = (0 - ct)t + 2a(7\ - T 2 ) 

•I ^MrT+S/rj (1 " exp ( ' m) - (35) 
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Now requiring that 

*,(T,) - tt (lib) 

results in 

•K - (0 - a)T 3 + 2a(7\ - T 2 ) 

(1 - exp (-T 3 /T,))(l - exp j-TJT,)) , , 
1- expt-dT. + TJ/r,] V ; 

which after some manipulation becomes 

r - (A - «)T 3 + 4a(7\ - T 2 ) / I coth ^ + coth ^r ' (37) 

Next, <f> 2 (t) is obtained from equation (35) by adding the term in 
equation (10) which corresponds to the transition time Cu = T*: 

2 (O - (Q - a)* + 2a(T 1 - T 2 ) 

T-ex P [-(r3 + r 4 )/r,] (1 cxp( ' /7l)) 

+ 2a(Z - T.) - 2a(T 1 - T 2 )(l - exp [-(/ - r 3 )/T,]). (38) 
Requiring that 

2 (T 3 + T 4 ) = 2ir (He) 

and performing some simple manipulation leads to the result 

2tt = Sl(T 3 + T 4 ) - a(T 3 - T,). (39) 

Equations (37) and (39) can be normalized by letting 

7 = a/a, (2) 

r,- = aT u i= 1, 2, 3, 4. (14) 

Equation (39) becomes 

2ir + t 3 — t 4 



7 ■ 



(15) 



r 3 + t., 

which is equation (15) in Section III. Equation (37) becomes 

y « 1 + i /, - 4<t, - r,)/ [coth IJ- + coth |j- }• (40) 

Eliminating y between equations (15) and (40) gives the constraint 
equation (16). 
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APPENDIX B 

The Phase Jitter Problem 

In this appendix it is proved that whenever the second-order PLL 
studied here achieves frequency lock, it also achieves phase lock: 

Lim <f>(t) = 2mr, | 7 I < Tp • (31) 

It can be shown that for this system any steady-state phase jitter #(£) 
is confined within the ±7r neighborhood of some lock point <f> = 2nir. 
This is a result of the periodicity of the phase-plane geometry in the <f> 
direction. However, since this proof requires a rather lengthy description 
of the properties of the phase-plane trajectories, it will be omitted. 
Below we prove that equation (31) holds when the phase error remains 
within such a ±ir neighborhood of a lock point. 

The technique used previously to calculate y p can also be employed 
here. Assume that in the steady state, a periodic phase jitter 4>{t) 
exists. Then the waveform f[<t>(t)] is binary and periodic, and the phase 
error is again given by equation (8). Since we are now considering a 
phase jitter within ±7r of <j> = 2mr, the requirements on <f>..(t) are: 

*..(0) = 4>..(T 3 ) - <t>..(T 3 + T 4 ) = 2717T. (41) 

Proceeding as in Appendix A, we obtain the equations 



2(r, - r 2 )(r a + r 4 ) = r 3 r 4 ^COth^- 4- coth ^- J • 



(43) 



These equations may be written directly from equations (15) and (16) 
by replacing ■* with 0. Below it is demonstrated that t 3 = t 4 = is the 
only solution of these equations when |y| < 1, which proves equation 

(31). 
Using |y| < 1 in equation (42) yields 

t 3 ,t<> 0, (44) 

so that only equation (43) must be considered. Using r of equation 
(18) and dividing by t 3 t 4 gives 



2r,(l - >•)(- +-) = COth;^- + COth^- 



(45) 
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Defining 

*t = V- , J = 3, 4, (46) 

equation (45) becomes 

(1 - r)(— + — ) - coth x 3 + coth x 4 . (47) 

Vr-3 X4/ 

Recalling that r ^ 0, we have for a; > 0: 

cc-thcO^^ 1 ^- (48) 

x a: 

Thus the only possible nonnegative solution of equation (47) is 

x 3 = x, = (49) 

which is the desired result. 
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